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In this paper, a new measure of entanglement for general pure bipartite states of two 
qutrits is formulated. 
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One main goal of modern quantum theory is the characterization and quantification of the 

in 
o 

key resource in many of the recent quantum information applications. In particular, quantum 
dense coding pQ, quantum teleportation [2j, and certain types of quantum cryptographic key 
distributions [Hill], rest crucially on the existence of entangled states. The question: "Given a 
quantum state, is it separable or entangled?", is nowadays completely solved for the simplest 
case of a pair of twodevel systems or qubits (i.e., 2x2 dimensional systems) [3]. Moreover, for 
this case, the Peres-Horodecki criterion of positivity of the partial transposition (PPT) [7j 
allows one to ascertain whether a general state of two qubits is separable or entangled. 1 However, 
concerning bipartite systems, no operational criteria have been established that allow to classify 
^ . a given state as separable or entangled for higher dimensional systems. On the other hand, 

several entanglement measures which quantify the amount of entanglement contained in a given 
state have been introduced in order to answer the following question: "Given an entangled 
state, how much is it entangled?" (for a review of the main entanglement measures and other 
related topics see, for instance, Refs. (SUniES^D)- For the special case of bipartite pure state 
PAB = IV'KV'I; a convenient measure of the degree of entanglement is the von Neumann entropy 
of either of the two subsystems A or B ^5] : 

E iW)) = -Tr(/?A log 2 Pa) = -Tr(p B log 2 p B ), (1) 

where pa (pb) is the partial trace of oy er subsystem B (A). The measure E(\t/j)) is also 

known as the entanglement of formation (EOF) of It can be shown that for the case of 
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1 The Peres-Horodecki criterion also applies to composite systems of dimension 2x3. For 2x2 and 2x3 
systems, the PPT is a necessary and sufficient condition for separability 0. 
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two qubits in the general pure state 
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1^2x2) = Yl 




(2) 



with normalization \aiij 



1, the EOF (1) is given by 0113 



£(|V>2x2» = h 
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(3) 



where h is the binary entropy function 



h(x) 



xlog 2 x — (1 — x) log 2 (l — x) 



(4) 



and where the concurrence C is 



C{\i>2x2)) = 2|ana 2 2 - 012021 



(5) 



E is a monotonically increasing function of C, and then the quantity C itself can be taken as a 
measure of entanglement. We also note that the concurrence (5) can be written equivalently as 



where k\ and k 2 are the two coefficients appearing in the Schmidt decomposition 15 , ^2x2) = 
^ll^ijyi) + K 2\X2-,y2) 1 and {|xi),|x 2 )} and {|yi),|?/2)} are orthonormal bases for the Hilbert 
spaces of subsytems A and B, respectively. 

In this paper, we formulate a measure of entanglement for a pair of three-level systems 
or qutrits (i.e., 3x3 dimensional systems) in an arbitrary pure state. This measure is the 
generalization to two qutrits of the concurrence for two qubits defined in either (5) or (6). Our 
measure of entanglement fulfills the following properties: (i) it is unique for a given quantum 
state; (ii) it is invariant under local unitary operations; (iii) it ranges in the interval [0, 1], the 
value (1) attained by the product (respectively, maximally entangled) state of two qutrits. 
Several aspects of qutrit entanglement have been previously considered by Caves and Milburn 
in Ref. JH], where the separability of various mixed states of two or more qutrits is investigated. 
Furthermore, in Ref. the degree of entanglement for a class of pure states of two qutrits is 
suggested in the context of a study of the Bell-CHSH inequality for two qutrits ^Hl IH3 EUl ■ We 
will then compare our measure of entanglement with the one proposed in ^2]; and show that 
the present measure constitutes the proper one characterizing the degree of entanglement for a 
general pure state of two qutrits. We mention that our work also represents an extension to two 
qutrits of the degree of entanglement for two qubits formulated in Ref. [U ■ 

We start by noting that any state pab of two qutrits can be expanded uniquely as ^S] 



C(\ip2x2)) = 2kik 2 



(6) 



PAB = 7: [ 1 ® 1 + \/3A A • 




(7) 



2 



where Aj (i = 1,2, ...,8) are the eight Hermitian generators of SU(3). For completeness, in 
the Appendix we write down the matrix representations of the generators Aj in the standard, 
orthonormal basis {|1), |2), |3}} for a qutrit. Both u = {m, . . . ,us} and v = {v\, . . . ,vg} can be 
regarded as vectors in a real, eight-dimensional vector space. Likewise, X A and X B are operator 
vectors acting in the Hilbert spaces of subsystems A and B, respectively. The (real) expansion 
coefficients in Eq. (7) are given by 

AT 

Ui = ^Tr(p AB \i®l), (8) 



■Tr(p AB l®\ j ), (9) 
Pij=^Tr(pABXi®Xj)- (10) 



2 
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It is not difficult to show that if pab corresponds to a pure state pab = IV") then we have 
that |u| = |v|. Furthermore, if \ip) is a product state then |u| = |v| = 1. This latter follows 
from the fact that the reduced density matrix of, say, subsystem A is given by 



Ti BPAB = l(l + V3\ A -u). (11) 



PA 3 

But, if pab = IV'KV'I an d \ip) is a product state, it is obvious that p\ = pa- Applying this 
condition to the operator (11), we obtain that |u| = 1. Note, on the other hand, that the case 
|u| = corresponds to the maximally mixed reduced density operators pa = Pb = gl- This in 
turn implies that the original pure state pab = IV'KV'I is maximally entangled. 
Let us now consider the scalar quantity 



c = VT^W. (12) 

The quantity C seems to be a good candidate to measure the relative amount of entanglement 
contained in a pure state pab = Indeed, as we have just seen, C = for the product state 

whereas C = 1 for the maximally entangled state. Admittedly, the main motivation to choose C 
as a measure of the entanglement present in a pure state of two qutrits is that the concurrence 
for two qubits in Eq. (5) or (6) can be equally defined by C(|V>2x2)) = \/l — |u| 2 , where in this 
case u is the Bloch vector determining the reduced density matrix pa = \ (1 + o A ■ u) |21| . 
Accordingly, we take the expression in Eq. (12) as the definition of the concurrence for a pure 
state of two qutrits, where the vector u is in turn given by Eq. (8). Thus, it can be shown that 
for a general pure state of two qutrits of the form 

3 

|^3x3> = ^2 (13) 

with normalization \dij\ 2 = 1, the concurrence (12) is given by 

CQipzxs)) = [3(|ana22 - ai 2 a 2 i\ 2 + |ai3«2i - ana 2 3\ 2 + |«i2a23 - ai3«22| 2 

+ |ai2031 - «lia32| 2 + |«11«33 - «13«3l| 2 + |a21<232 - 2 2a3l| 2 

+ |«23031 - «21«33| 2 + |«13032 ~ ai2«33| 2 + |a22033 ~ a23032| 2 )] 2 • (14) 
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Eq. (14) is the generalization to two qutrits of the concurrence for two qubits in (5). It is helpful 
for what follows to introduce the matrix with entries given by aij in (13), 



a — 0*21 «22 «23 • (15) 




It is a nice elementary exercise to show that if 1^3x3) is a product state then necessarily det a = 
(please note, however, that the converse of this statement is not true). In the same way, it is 
readily seen that if \1lj3x3) is a product state then each of the nine terms appearing in Eq. (14) 
vanishes, so that C(|V>3x3)) = 0. On the other hand, if \1lj3x3) is a maximally entangled state 
then the matrix a is restricted to obey the conditions pa = = |l and pb = a^a = gl. 
As a result, the coefficients a^'s are constrained in such a way that they cause C(\ip3 X 3)) 
to attain its maximum value C(\ipsx3)) = 1- We note, incidentally, that from the condition 
&&* = |l or a^a = |l one obtains immediately the result that, for the maximally entangled 
state, I deta| 2 = 27- Of course, for an arbitrary state the concurrence (14) lies in the interval 
< C(k/>3X3)) < !• 

The state in (13) can be written alternatively in terms of a Schmidt decomposition ^3] 

1^3x3) = «i|a?i,2/i) + K 2 \x 2 ,y2) + K 3 |x 3 ,y 3 ), (16) 

where and are two orthonormal basis sets belonging to the Hilbert spaces of sub- 

systems A and B, respectively. k±, K2, and re 3 are real and nonnegative coefficients satisfying 
k| + k| + 4 = 1 (without loss of generality we may also take K\ > K2 > ^3)- Since k±, k 2 , and 
k 3 are unique for any given state \tp3x3), it should be possible to uniquely express C(|?/>3x3)) i n 
terms of the Schmidt coefficients K\, k>2, and k 3 for the state 1^3x3)- From the theory of the 
Schmidt decomposition ^S], we know that the squares k\, 4i an d 4 correspond to the three 
eigenvalues of the reduced matrix pa = aa^ or ps = aVa. These eigenvalues are solutions of 
the cubic equation 

A 3 - A 2 + Ic 2 (|V> 3x3 »A - I det a\ 2 = 0, (17) 

with C(|V>3x3)) being the expression in Eq. (14). The roots of the equation in (17) are real and 
nonnegative. Denoting such roots by At 2 , ac 2 , and k|, it follows from the coefficients of the cubic 
equation (17) that its roots satisfy the relations: 

k\ + k\ + k\ = 1, (18) 
44 + 44 + 44 = ^C 2 (\^ 3x3 )), (19) 
K i K 2 K 3 = I det a| 2 . (20) 

From Eq. (19), we obtain the concurrence C(\^3 X 3)) m terms of the Schmidt coefficients asso- 
ciated with 1^3x3) 

(21) 



c(|^ 3 x3» = ^3(44 + 44 + 44). 
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Eq. (21) is the generalization to two qutrits of the concurrence for two qubits in (6). As we have 
seen, for the maximally entangled state we have C(|V>3x3)) = 1 and | deta| 2 = Substituting 
these values into Eqs. (18)-(20) gives k\ = k| = k\ = |. On the other hand, for product states 
we have C(\ip3 X 3)) = and deta = 0, which gives k\ = 1, k\ = k| = 0. The case in which 
C(|V>3x3)) 7^ and deta = 0, corresponds to the states with k\,k^ ^ and k\ = 0. When one 
of the Schmidt coefficients (say k 3 ) is zero, the concurrence (21) reduces to C(|V'3x3(^3 = 0))) 
= ^/3kik 2 , ranging in the interval [0, ^\/3]. Therefore, the EOF for the class of states 

1^3x3(^3 = 0)) is 

S(lfe3(*3 = 0))) = A ( 1 + ^™ ). (22) 

£"(1^3x3(^3 = 0))) ranges from for the product state (i.e., C = 0) to 1 ebit for the "singlet" state 
\1p3x3} = (l/v / 2)(|2;i, Hi) + 1 3^2 , 2/2 ) ) (i-e., C = \/3/2). Naturally, we know that the maximum 
EOF is obtained for the maximally entangled state 1^3x3) = (1/ \/3) (|^i , J/i ) + 1^2, J/2)) + |a?3 , 2/3) ) ? 
this maximum value being equal to log 2 3 ~ 1.58 ebits. Unfortunately, however, we lack a general 
formula for the entanglement of formation -E (1^3x3)) of an arbitrary pure state 1^3x3) m terms 
of the concurrence (21). 

Of course, the entanglement measure (21) satisfies the property of invariance under local 
unitary operations. This follows from the facts: (i) and the eigenvalues of either 

the matrix pa = &6r or ps = a^a; and (ii) The eigenvalues of a matrix are unaltered if the 
matrix is transformed by a similarity transformation. In particular, this latter applies to local 
unitary transformations UaPaU a and UbPbU b , where Ua and Ub stand for general 3x3 unitary 
matrices acting in the Hilbert spaces of subsystems A and B, respectively. 

We conclude by recalling the degree of entanglement for two qutrits proposed by Fu et al. 
[TTj . and then comparing their measure with ours. Fu et al. considered the class of states of two 
qutrits in the (Schmidt) form 

= -J=(ai|l,l)+a 2 |2,2) +a 3 |3,3)), (23) 

where the a^'s are real coefficients varying in [— \/3, v3], and satisfying a\ + a\ + a| = 3. By 
studying the Bell-CHSH inequality for two qutrits (with two von Neumann measurements per 
site) [H3 033 120] j S < 2, where the specific form of the quantity S is given in Eq. (5) of [T7| . 
Fu et al. showed that, for the class of states (23), the quantum-mechanical prediction of S 
acquires the form 

SqmG^)) = aia 2 T 12 + 0103^13 + a2a 3 T 23 , (24) 

where the Tjj's depends solely on the local measurement settings. Furthermore, Fu et al. proved 
that the maximum value of S'qmGV')) is determined by the absolute value of the cross products 
a\a 2 , 0103, and 0203. In view of this, Fu et al. suggested to describe the degree of entanglement 
for the two-qutrit state (23) by 

Mm = M + M + M ^ (25) 
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Fig. 1. Degree of entanglement of the state (23) for the case in which a\ — 0,2, where 
Pe and C are defined by Eqs. (25) and (26), respectively. For a% — 02, the state (23) is 
completely characterized by a single parameter e. 

On the other hand, for the state (23), the concurrence (14) reads 

C(W)=( |QlQ2|2 + |Ql ^ |2 + |Q2Q3|2 y. (26) 

In Fig. 1, the functions Pe(\iP}) and C(|V>)) have been plotted for the case in which a\ = ai = 
a/36/2 and a 3 = ^3(1 - e), where < e < 1. We can see that C > Pe for the entire range of 
variation of e. For e = (i.e., for the product state) we have Pe = C = 0, whereas for e = | 
(i.e., for the maximally entangled state) we have Pe = C = 1. The case e = 1 corresponds 
to the singlet state \ip s ) = (l/v / 2)(|l, 1) + |2, 2) ) , and for this case Pe = \ and C = |\/3. It 
is to be noted that the value Pe = \ is too low to adequately describe the relative amount of 
entanglement present in the state |^ s ), as such value leads to an underestimated entanglement of 
formation of the singlet state. Indeed, substituting C by \ in Eq. (22) would give E(\ip s )) 0.44 
ebits. However, we know that the EOF of \ip s ) is, by definition, equal to 1 ebit ^S]- We thus 
deduce that the correct measure of entanglement over the entire range < e < 1 (and, in fact, 
for any state (23)) is the function in Eq. (26). 

To finish, we consider the generalization of the concurrence in Eq. (21) to the case of a pair 
of d-dimensional systems or qudits {d > 2) in a pure state. Let us assume that the state of the 
two qudits is written in the Schmidt form 

\^dxd) = Ki\xi,yi) + K 2 \x 2 ,y2) H h K d \x d ,y d ), (27) 

where and {|yj)} are orthonormal bases for the Hilbert spaces of subsystems A and B, 

and the «j's are real and nonnegative coefficients satisfying k\ + k\ + • • • + = 1. Then the 
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concurrence for the state (27) could be defined by 




C(|r f/X ,,;,) - s : j -f^j ( K 2 K 2 + ... + K 2 K 2 + K 2 K 2 + ... + K 2 K 2 + ... + K 2_ iK 2) { 2sS ) 



Note that expression (28) reduces to the concurrences (6) and (21) for the cases d = 2 and 
d = 3, respectively. Moreover, C(|^xd}) = f° r the product state, whereas CQijjdxd}) = 1 for 
the maximally entangled state, i.e. the state in (27) for which K{ = \ j\fd~ for each i = 1, . . . , d. 

In summary, in this paper we have presented a new measure to quantify the degree of 
entanglement for two qutrits in a general pure state |^>3x3)- This measure has been expressed in 
two equivalent forms: either in terms of the expansion coefficients ctij of \1p3x3) in an arbitrary 
basis \i,j), i,j = 1,2,3 (cf. Eq. (14)); or else in terms of the coefficients Ki of |^3x3) in the 
Schmidt basis {\x{, yi}}, i = 1, 2, 3 (cf. Eq. (21)). We have derived the entanglement of formation 
of 1^3x3) for the simple case in which one of the Schmidt coefficients is zero. Also we have 
compared our measure of entanglement for two qutrits with the one previously proposed in 
Ref- |17j . Finally, we have extended the concurrence for two qutrits (21) to the case of two 
qudits described by the pure state (27). 



Appendix 

In the orthonormal basis {|1), |2), |3)}, the generators of SU(3) have the matrix representations: 
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The generators Aj obey characteristic commutation and anticommutation relations (see, for 
example, Refs. |16l 1221 123j). Here we merely note that the matrices (Al) are traceless and 
satisfy 

Tr(X a \p)=25 aP , a,/3 = 0,...,8, (A2) 



where we have supplemented the eight generators Aj with the operator Ao = y 1 1 |16j . Relation 
(A2) is useful in deriving the coefficients (8)- (10). 
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